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During : the six months from April to October, 1985, progress was made 
with both the Euler code calculation task and the rotor optimization task as 
briefly described below, 

EULER CODE TASK 

During this period a number of refinements in the computer code were 
worked out and tested. Three codes have been written to date. One program 
is for am isolated wing and is being used to compare with data for the vortex 
wake (Weston) . The second code is for an isolated wing with a streamwise 
vortex passing above it. This program is being used to validate the com- 
putational procedure for incorporating the vortex into the Euler equation 
calculations. The third program is the hovering rotor code which is the 
overall objective of the research. In this six month period, attention has 
been given mostly to the first two of these three programs. An attached 
abstract submitted to the AIAA 4th Applied Aerodynamics Conference in San 
Diego in June describes our progress on these programs. 

We plan to turn our attention to' the third program to get it updated 
during January 1986, and then to perform the final calculations of the study 
during February 1986, A more detailed description of our progress is being 
written and will be available during the next period. 

ROTOR OPTIMIZATION TASK 

The optimization calculations for a hovering helicopter rotor have 
been completed and are being prepared in a thesis by Song Young Chung 
which will be submitted later. 
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Euler Solutions for the Flow Around a 
Hovering Helicopter Rotor 
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1 Introduction 

In recent years, computational fluid dynamics technology has been highly 
developed for fixed wing aircraft configurations. The same cannot be said for 
rotary wing aircraft. In the case of a helicopter rotor in hovering flight, which 
may be treated as steady in the frame of reference of the rotor, an important 
effect is the interaction of the tip vortex trailed from one blade with the following 
blade. This interaction results in large variations in the loading at the blade 
tip. Furthermore, the semi-infinite vortex wake of the rotor has a significant 
influence on performance and must be included. The need to include these 
effects poses quite a challenge for existing computational techniques. 

The computation of aerodynamic flows using the compressible Euler equa- 
tions allows the “capturing” of vortical wakes trailed from lifting surfaces. How- 
ever, due to grid stretching and numerical viscosity, the wake vorticity is diffused 
as it is convected downstream. For this reason, an Euler calculation around a 
complete rotor configuration is impractical, as the tip vortex of each rotor blade 
will be quite diffuse by the time it reaches the following blade. The steep span- 
wise gradients in the aerodynamic loading near the tip due to the blade/vortex 
interaction will not be correctly computed. Also, a model for that part of the 
semi-infinite wake that lies outside the computational domain is required if ac- 
curate values for the thrust, induced power, and the initial radial contraction 
of the wake are to be obtained. 

In the proposed paper, we will solve for the flow field around a hovering 
helicopter rotor using an approach that contains the following features: 

1. an Euler calculation of the near field of a rotor blade, including the for- 
mation and roll-up of the attached vortex wake; 

2. treatment of the tip vortex of the proceeding blade such that the vortex 
is not destroyed by numerical diffusion; and 

3. inclusion of the influence of the semi- infinite wake. 

2 Technical Approach 

The basic scheme for solution of the Euler equations for the flow around 
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the rotor blade is a finite volume, multi-stage algorithm based on the work of 
Jameson and Baker [l). An 0-0 mesh topology is used to get good resolution of 
the tip, leading edge, and trailing edge regions of the blade (fig. 1). The semi- 
infinite vortex wake is treated using by the free wake model of Miller [2,3]. The 
coupling of the free wake and the near field blade calculations will be performed 
in a fashion similar to the approach used by Roberts and Murman [3] for a 
lifting line representation of the blade, and is outlined below. 

To account for the close encounter between the rotor blade and the tip vortex 
from the proceeding blade a prescribed vortex, or perturbation, approach will be 
used to include this vortex without the attendant numerical diffusion problems. 
This perturbation approach has been applied to the unsteady 2-D blade/vortex 
interaction problem by Srinivasan and his co-workers [4j. The implementation of 
the approach hero consists of freesing the wake geometry, including the position 
of the tip vortex of the proceeding blade, during the Euler iterations. The state 
vector of the wake flow field, Uq = (p, pu, pv, pw, pE) is then computed. The 
finite volume flux balance and and artificial viscosity operator are applied to this 
flow field. The residuals of the flow field Uq are subtracted from the residuals 
of the total flow field U at each time step of the Euler iteration. In this way, the 
truncation error of the scheme is corrected for the influence of the tip vortex of 
the proceeding rotor blade such that the vortex remains well defined, even near 
the outer boundary of the computational domain where the resolution of the 
grid is quite coarse. 

With the new Euler solution, the strengths of the wake vortices are re- 
computed, and the wake geometry is updated through the free wake iteration 
procedure. After updating the wake geometry, the new geometry is used to re- 
compute the prescribed flow field Uq, and the Euler solution is updated. This 
coupled iteration procedure is continued until convergence. 


3 Results 

Euler solutions have been obtained for fixed wing geometries. A comparison 
of solutions obtained with this code and experimental data for the ONERA M6 
wing at a = 3.06°, Moo = .84 are shown in figure 2. The calculations were 
performed on a grid of 96 x 20 x 20 finite volume cells, and comparisons of 
surface pressure coefficient for 3 span stations are shown. 

The code has also been used to compute the flow about an unswept, unta- 
pered wing with a NACA 0012 airfoil section tested at NASA-Langley Research 
Center by Weston [Sj. The purpose of this calculation is to compare the vortical 
wakes computed by the Euler equations with experimental data. The calcula- 
tion was performed on a 128 x 32 x 32 grid at a = 8°, Moo = .1425. Surface 
pressures are shown in figure 3. In figure 4 are shown contours of constant total 
pressure and vorticity at 1/2 chord behind the trailing edge. The computed 
wake is qualitative similar to experimental data, indicating that the Euler equa- 
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tions do capture the vortical structure of the wake. Quantitative comparisons 
with experiment will be performed to ascertain the ability of the Euler equations 
to capture realistic vortex wake structures. 

The use of the perturbation approach for computing the interaction of a 
stream wise vortex with a semispan wing is shown in figures 5 and 6. The ge- 
ometry of the configuration is sketched in figure 5a; the wing has a semispan 
to chord ratio of 2:1, a NACA 0006 airfoil section, and the free stream Mach 
number is .3. The vortex passes 1/2 chord above the wing and 1/2 semispan, 
from the plane of symmetry. A 96 x 20 X 20 grid was used around the wing. The 
spanwise distribution of lift coefficient Ci is shown in figure 5a for the perturba- 
tion approach, while figure 5b shows the results for the standard Euler solver. 
The steeper gradients at the 1/2 semispan section for the perturbation approach 
due to the better resolution of the vortex are apparent. In figure 6, contours 
of vorticity are shown in a vertical plane passing through the mid-chord of the 
wing. The contour levels are identical in figures 6a and b. Note that the vortex 
is well defined with the perturbation approach, while for the non-perturbation 
scheme it has virtually vanished. This calculation illustrates the ability of the 
perturbation scheme to overcome the problems of coarse grid resolution near 
the outer boundary and the diffusion of vorticity due to artificial viscosity. This 
is a preliminary calculation; further comparisons with the experimental data of 
reference [6] for this configuration will be undertaken. 

For the proposed paper, the current perturbation scheme will be applied to 
the case of a hovering helicopter rotor. The Euler code for a hovering rotor has 
been written. Preliminary results for the hover case have been presented by 
the authors [7] using an earlier version of the perturbation scheme now being 
developed. (In reference [7], only the smoothing operator and not the flux 
operator was applied to the wake flow field Uq.) The position of the tip vortex 
passing near the rotor blade will be found from a free wake calculation. The free 
wake and Euler computations will be coupled such that the strengths of the free 
wake vortices will be determined from the Euler calculation. In this fashion, 
the correct representation of the tip vortex of the proceeding blade within the 
Euler computational domain as well as the influence of the entire semi-infinite 
wake will be included in the calculations. 


4 Summary 

A proscribed vortex or perturbation scheme for the Euler equations around 
a hovering helicopter rotor is being developed. For the rotor case, this approach 
will provide several advantages over existing hover performance prediction meth- 
ods. First, the scheme will allow the interaction between the rotor blade and the 
tip vortex of the proceeding blade to be computed. Second, the formation and 
roll-up of the near trailing vortex wake attached to the blade will be captured 
by the Euler solver. Finally, by fully coupling the Euler solution of the near 
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field of the blade with a free wake analysis, the proper influence of the complete 
semi-infinite wake will be accounted for. 
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Figure 1 - 0-0 Mesh Features 
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Figure 2a - C p Comparisons for ONERA M6 Wing - 20% span 


ONERA M6 — 3.06 alpha, .84 Mach, 96 x 20 x 20 grid 
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Figure 2b - C p Comparisons for ONERA M6 Wing - 65% span 
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M6 Wing - 90% span 
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WESTON TEST CASE — 8 alpha, .1425 Mach, 128 x 32 x 32 grid 
8 CHORDWISE SECTION 



Figure 3c - C p for Plain, Unswept Wing - 
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Figure 3d - C p for Plain, Unswept Wing - 99% span 
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